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The theory of plasma waves and Landau damping in Maxwellian plasmas, Landau's "rule of pass 
around poles" include doubtful statements, particularly related to an artificial "constructing" of the 
dispersion equation, what should allow the possibility of its solution otherwise not existing at all, 
and the possibility of analytical continuations of corresponding very specific ruptured functions in 
the one- dimensional Laplace transformation, used by Landau, what is the base of his theory. 

We represent, as an accessible variant, a more general alternative theory based on a two- 
dimensional Laplace transformation, leading to an asymptotical in time and space solution as a 
complicated superposition of coupled damping and non-damping plane waves and oscillations with 
different dispersion laws for every constituent mode. This theory naturally and very simply explains 
paradoxes of the phenomenon of plasma echo. We propose for discussion a new ideology of plasma 
waves (both electron and ion-acoustic waves) qualitatively different from the traditional theory of 
Landau damping for non-collisional as well as for low-collisional plasmas. 



The sophisticated theory of Landau damping of longitudinal plasma waves, which was published by Landau in 
1945 [1], is now considered as the greatest success of the physics of plasma, the heritage verified experimentally and 
theoretically, which has not to be doubted. But this halo is to a considerable extent dispersed after some critical 
analysis of different sides of the problem. 

1. According to the classical linear theory of longitudinal plasma waves, in order to derive the dispersion equation 
one has to substitute expressions of the type of exp(— iuit + ikx) for electrical field E and perturbation fi(v, x, t) of 
the isotropic electron distribution function fo(v) into the linearized non-collisional kinetic equation 



and Poisson's equation 



^Il +Vx ^Il + e Jk^Il = Q (1) 

dt x dx m e dv T ' 



where E x is the selfconsistent electrical field in x direction, v x is the ^-component of the electron velocity, e, m e 
are correspondingly the electron charge and mass, u> and k are generally complex; and then to solve the resulting 
dispersion equation D(tu, k) — defining a dependence k — f(u>) for electron plasma waves, one or more, according 
to the number of its roots. 

For plasma waves, as it can be immediately shown assuming oj = luq — id at real u>q, S and fc, breaking up integrals 
according to velocity components in the range — oo < ?J, < oo to constituents — oo < v x < and < v x < oo and 
picking out the real and imaginary parts, after some very simple procedures: the dispersion equation has no solutions 
for S. So, one obtains easily for D(lu, k) and imaginary part Im_D(oj, fc) = 0, determining preferably (5, equations 

nf n i 47rze 2 f df (v) dv _ 

D{u,k) = l-— — / — — = l-iF(p), 

km J_ 00 ov x p + ikv x 
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lmD(UJ > k) = L dVy L dVz L d^H-vvi-^+As^r (3) 

where D(uj, k) is susceptibility, F{p) is a particular case of "Landau function" , which will be discussed later, p = —iuj, 
and dfo/dv x < 0. Since near the point v x ~ coo/k the integral appears to diverge as 1/(5, the whole expression 
for Imi)(w,fc) seems to be finite, though ruptured at S — > ±0. But here the integrand in the left-hand side is sign- 
invariant for both Maxwellian /o and any other /o with sign-invariatn dfo/dv x , so this equation is evidently impossible 
to satisfy at whatever iv ^ and any 5, including S — > ±0 (in contrary to ReD(w,fc) = 0). Strongly in 5 = this 
expression makes no sense. This equation has no solutions, too, independently on Landau theory and his analytical 
continuations of such type of functions, on sophisticated Van Kampen modes and whatever theories else. 

The case of a low-collisonal plasma is analogous, but with replacing S, in the simplest approximation, by (5 — v), 
where v is collision frequency, with the same result of non-existing solutions. 

This statement is also equivalent to the note in [1] about non-existence of complex poles in the right halfplane of the 
Laplace transformation used in [1]. In the very strong canonical sense it means indeed non-existence of solutions of 
the type of a single (or of some number, according to possible roots at fixed k, as in multimode Landau's asymptotical 
solution) monochromatic plasma waves. 

One may also expand the terms of the dispersion equation, after integration over velocity spherical angles dtpdO, in 
the small parameter 5/ui n <C 1, but the resulting series is asymptotically divergent (e.g., cf. [2]), and according to the 
theory of such expansions, the error of a resulting (in this approximation) dispersion equation solution for S (which 
just coincides here with the traditional Landau's expression for S) turns out to be comparable in value with the error, 
defined by a linear in S term of the expansion in 5, i.e. the solution error is comparable with this solution itself. 

2. Landau's the most slowly damping asymptotical wave should indeed satisfy the dispersion equation (3). But the 
existence of a solution of the disperation equation is attained by including in it some additional terms defined according 
to considerations not related to the dispersion equation insclf. It is found that such "artificially constructed" (i.e. 
to the Vlasov integral in the arbitrary sense of principal value one arbitrarily adds the half-residium with semi-circle 
contour in the complex u x -plane in the sense of Landau's rule) dispersion equation which just now has solutions, 
corresponds at asymptotically large times to plane attenuating waves. But it is trivially evident that the dispersion 
relation does not contain variables x and t, and at t — > oo the precise dispersion equation does not turn into the 
modified dispersion equation of Landau theory. 

The terms added by Landau [1] are found by a calculation of some contour integral along the real axis in the plane 
of the complex variable v x with passing around the pole —iS/k near the point u a /k of the real axis at S — > +0 (Fig.l). 



-iS 



Fig. 1 (a) and (b). Some alternative variants of calculation (with different results) of the Landau contour 
integrals F(p), p = —iu, in the plane of complex v with u> = wq — iS (real u>o, 5 > 0) at S — > +0. 



It is easy to see that in dependence on a way of tending S to zero there are equally accessible, as an example, 
both possibilities: (a) (half-residium in a pole near the real axis, what is a result of the Landau rule of going around 
the pole), and (b) (total residium in the pole), as well as other variants, including those leading to possible infinite 
results. Any discussion about a preferable Tightness of either the first or the second result (see e.g. [3]) appears to 
be objectless because of non-existing any definite limit of the contour integral at 5 — > at all. This means, that the 
result of this calculations depends on a selected way of tending 5 — > relative to the contour integral calculation 
procedure in the complex plane v x . Because of it an analytical continuation of the regular in the upper halfplane of 
complex lo = ip function F(p) into the lower halfplane w appears to be non existing. 

For calculation of expressions of the type (3) at 6 — ► ±0 one might use the relation which is known in mathematics 
as Sokhotsky formula 



lim 

e^±0 



f°° f(x)dx rfn . f°° f(x)dx 

/ ^ = T«/(0) + v.p. (4) 



which thus appears to be only a particular case because of the limit existence being related to a selected definite 
procedure of calculation of an integral (according to the definition) as a limit of the corresponding sum 
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1-°° fjx)dx = Um yji^l^ 

J_ oc x + ie Ax h ^a^x k + ie 

Substitution of this definition (5) into the left-hand side of eq. (4) results in a double-limit sum, the value of which, 
as being a limit, is evidently depending on the arbitrary choice of the relative cooperative rates of tending to zero of 
the values of Ax k and e. 

Strictly speaking, it takes plase also for the above considered expression (3) lim[Im£>(w, k)} (and F(p)) at S — > ±0. 
So, the above pointed out result of some definite finite limits for that one is only some particular way of calculation of 
the mentioned double-limit indefinite sum. That is, the expression \im[Im D(w, k)] in (3) has really no definite sense 
at all. 

3. In the well known paper of Landau [1] the field function E(x,t) was found by solving coupled equations (1) 
and (2) at an in advance arbitrarily taken coordinate dependence E{x,t) = <p(t) .exp(ikx) and /i ~ exp(ikx) at all 
the times (cf. also [4], [5]) by means of the one-sided one dimensional Laplace time transformation with the use of an 
analytical continuation of a function of the form 

from the right half-plane of complex p to the left through the peculiar line Rep = e = and with the plane damping 
wave as an asymptotical solution (v is hereinafter x-component of velocity with omitted subscript x; ip{v) is a real 
function of v). 

Introducing complex v with the counterclockwise rule of going around the pole v p = ip/k at Rep < (Landau's rule) 
at integration by v, one removes the rupture at the line e — > ±0 defined e.g., as some partial cases, by the Sokhonsky 
formula (half-residium) or by the residue theorem, but does not remove the rupture of derivatives dF(p)/dp in the 
analyticity Cauchy-Ricmann conditions in the direction which is perpendicular to the line Rep = 0. So, at Rep < 
one finds, if, e.g. using the residue theorem in the calculation of the contour integral (6), correspondingly at the 
left-hand side (— ) and at the right-hand side (+) from the line Rep = 0: 

(that is the Landau's type analytical continuation to the left), and 

*«-£££■ <■>"■ 

At the calculation of derivatives across (perpendicularly to) the line Rep = with a function tp having not 
peculiarities, one finds a finite expression 

= (e^-0), (9) 

but contrary to it a divergent (infinite) at e — > +0 integral expression (with infinite principal value at e = 0) for the 
derivative of the analytical at e ^ function (8) 



dF+(p) _ f°° ip(v)dv 



dp J_ 00 (p + ikv) 
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(e-+0). (10) 



4. To the very hard perceived consequences of the Landau theory one should attribute also the existence of 
the pole p = —ikv, which appears at the Landau procedure of analytical continuation of the Laplace transform of 
the distribution function fi(v,x,t). This pole leads to an additive term for the distribution function of the type 
H(v) exp(ikx — ikvt), which one uses for an interpretation of the plasma echo [4]. Dramatic attempts to comprehend 
and to conciliate this "mysterious" term with the common sense are clearly demonstrated in the textbook [4] and 
are related to the appearing paradoxicability of the existence of non-damping plasma oscillations in the absence of 
whatever recovery force, whereas the electric field had to be disappearing due to the Landau damping! 

Such /i, independent on E(x,t), can not satisfy the Poisson equation (2), since its substitution in the right hand 
side of (2) leads to quite arbitrary functions of time and is not consistent with the field damping in the left hand side. 
At the same time the paradox of plasma echo appears to be naturally explained in principle in the further developed 
theory. 
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5. In analogy with [1] let us consider a semiinfinite (halfspace slab) plasma with an initial perturbation 

/.<-Ho (K1,exp(,M ::<:, <"> 

where v is now the velocity component along x, and is some integrable function of \v\. 

Equation (1) is easily solved according to the known procedure with the method of characteristics (by reduction to 
the equivalent system of simple ordinary differential equations) with the result 

fi{v,x,t) = g{\v\)e ikx + [ w[x-v(t-t'),t']dt' (12) 

J (x,x'>Q) 

where 

eE(x,t)df . . 

wM ^-^r^ (13) 

x' = x - \v\{t - t') (14) 

and an additional physical condition x' > is related to the finite speed of spreading of particles which are arriving 
to the point x. 

Substituting expressions (12) - (14) into eq. (2) one obtains a linear integro-differential equation for E(x,t): 

/OO f-t 
dvv exp(-m e v 2 /2k B T) x / dt'E[x - v(t - t'), t'] - 47rea(x) exp(ikx) (15) 
-OO Jtn(v) 



dE(x, t) ^/8-Km e e n e 



dx k B T 



where the dependence e~xjp(ikx) is supposed only for the sake of correspondence with the assumed by Landau coordinate 
dependence [1], and a(x) = a = const at x > ; a(x) = at x < ; t n (v) = at t — x/\v\ < ; t (v) =t — x/\v\ at 
t-x/\v\>0. 

The one-value solution of eq. (15) must be defined in the same extent by both initial and boundary conditions. 
The dependence E(x,t) on t, as well as on x is defined by solving equation (15), and some analogy with the Landau 
problem [1] for the infinite plasma may be realized by taking the limit x — > oo. 

The general solution of equation (15) can be found with onesided two-dimensional Laplace transformation: 

E((,t) = -—2 e^dpj E pip2 e^dp 2 ; (16) 

yZTTl j J a\—ioo J o^— ioc 

8( (2m) 2 J ai _ loo J a2 - ioo 27TZ 7 ff2 _ JOO V2 

where o\ , a 2 > , and there are introduced normalized dimensionlcss variables 



' / 777 I 

C = kx > T = t/ta > ta= ki2kif = k^' (18) 

k is real; vt is thermal mean velocity. Here appears also some E(0, t) as a boundary field constant of integration. 
It can be naturally considered as a given field related to some external circumstances as it is due to some external 
sources, etc. 

After substitution of eqs. (16), (17) and (18) into eq. (15) one obtains 



-47r 2 ai(C)e lC +,4 / d Pl dp 2 E PlP2 e piT+p ^ x 

J <T\— iOQ J <72—iOQ 



fC/ T 2 I _ e ~p2^r-p 1 T l_ e P2£,r-p 1 r roc ^ ^ _ ^- pi Q / ^- p2 Q y — e ~PiC/C+P2C 

{ / d ^ I — — c Z~ 1 + / d ^ I — — c 7 ]} 

Pi+P2? Pi-P2t, Jc/t Pi+P2t, P\-piK 



dpi j dp 2 p 2 E pip2 e^ T+p ^ -2itiE(0,t) / — dp 2 , (19) 



o 

<T2— tOO P2 

where A = const; stepwise function ot\ is proportional to a(x). 
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Keeping in mind the subsequent consideration of some longitudinal analogy of the skin effect, we shall suppose the 
field E(0, t) takes the form 



E(0,t) =E o e- i0T , (20) 



dpi, (21) 



where (5 is some real number. 

Taking in account equalities of the type 

e -i0r _ 

2 ™ Jai-ioo Pi +iP 

and finiteness of the integrands in points £ with p\ ± p2£ —> for terms of the type 

I _ g±(j>l±P2£) 

Pi ± P2^ 

so one assuming the following integrations along £ in the sense of the principal value, one can find an asymptotical 
solution of eq. (19) as 

a" l -Eq 

_ PlP2(p2-i) p 2 2 (pi+i!3) . . 

where a* ~ ao , and £ is some characteristic value of the normalised velocity £ = v/vt; £ and A£ are of the order of 
1, and 

< 23 » 

In the derivation of eq. (22) it is supposed that one may apply to integrals over £ in the braces of eq. (19) the 
mean value theorem, replacing functions of £ in the integrands with the functions of some constant value £ belonging 
to the interval of integration. Also at integration over pi , p 2 one can neglect arising terms 

„-p 2 ?T-piT p -PlC/C-P2C 

— ■ (24) 

Pl + P2t, Pi + P2^ 

asymptotically for <j\, oi > and large £, r — > oo. At some fixed value of (72 > Laplace transformation implies the 
integration contour o\ ± zoo to be at right of the pole p\ — £>2£ = (<ti > CT2£); this also means the possibility to 
neglect in an asymptotical solution with large £, t — > oo the exponentially small terms 

pP2^T-p 1 T p -PlC/C+P2C 

j\ ~ J- (25) 

P1-P2K P1-P2K 

Asymptotical value of E((, r) at both ( — > 00 and r — > 00 is then defined by the pair poles of expression (22) at 
shifting the integration contours (71,02 to the left. To the term with ao ^ there corresponds a constituent of the 
asymptotical solution £?(£, r), defined by the residua sum in the poles 



(Pi = 0, p 2 = i); (pi = 0, p 2 = ±^2Ae-i 2 A0; 



(pi = i^Ae-^AC + l, p 2 = *) (26) 

(the poles (pi,fJ2) = (0, 0) do not contribute to E((, t) since E PlP2 — > at successive transits p\ — > 0, then ^2 — > or 
v.v.). 

The phase speed of waves, which correspond to the last pair (pi,P2) in (26) at small n e , coincides with the mean 
thermal velocity 



wo _ J_ = / 2fc B r 

fc fofc V 
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what is in accordance with experimental plasma echo speeds. At large A one finds 



wo ^ V^A _ ^ wi 

where w^, if Langmuir frequency, and 7 is a factor of the order of 1 . 

For the term in (22) with Eq 7^ (with expected extinguishing plane longitudinal waves) one obtains pole pairs 
defining additional constituents of the asymptotical solution: 



(pi = -i/3, P2 = 0); (pi = ±iM2Ae-t 2 A£, p 2 = 0); 



(pi - -i/3, P2 = ±^/2^le-« 2 A^ - (/3/0 2 ). (27) 

But in the pole p 2 = of the second order (that is l/p\) the residium equals zero. 

The "mixed" constituent of solution which corresponds to the poles (pi = —i(3,p2 = i) at x,t — > 00 , gives zero 
contribution. 

The derived asymptotical solutions content some exponentially growing modes. This is easy explained by the non- 
selfconsistency of assumed initial and boundary conditions fi(v, x,0) and E(Q, t) in the absence of an outside source 
field. The divergent terms can be easily removed with the corresponding fitting amplitude values a* and Eq in eq. 
(22) and assuming /3 = 0. So one obtains only nondamping standing oscillations as a sum of modes moving with the 
Langmuir velocity at large n e and the one-dimensional mean thermal velocity vt at small n e . 

One may suppose that an outside source field might be accounted for by dividing the boundary value of the field 
into a part of the selfconsistent plasma response and the rest one. The first one is defined by condition of the mutual 
cancellation of terms, which exponentially grow at x — ► 00. The second one gives proper plasma modes induced by 
an outside field. 

The above found very simple result of a non-damping standing wave for above variant of the self-consistent initial 
and boundary conditions seems to be very natural physically. 

In a variant of a pulse exitation of longitudinal waves with the field E along x at 

a(x)=0; E(0,t) = E 5+(t) (28) 

one can obtain with the Laplace transformation of S- function 5+(t) being 1 (see [6]) a solution 

E - E ° /l>2 (29) 



-, 2Aee-Z 2 AZ ' 

1+ Pi-vie 



respectively only two pair poles 



(pi = i^2Ae-t 2 AZ, P2 = 0); ( Pl = -i^2Ae~i 2 A^ p 2 = 0); (30) 

with zero residua and asymptotic value E(x,t). 

These results, which are in principle accessible to experimental verifications, are discrepant qualitatively with the 
Landau theory, so this means the necessity of a new ideology of plasma waves. The asymptotic limit is represented 
by a superposition of a finite number of coupled both damping and non-damping oscillatory modes with different 
dispersion laws. 

Using the method of Laplace transformation one may easily obtain also asymptotic solution for the perturbation 
fi(v, x,Q). In this case an arising pole (p\ + P2O and others may indeed lead among others to terms in /1 of the 
type exp[ik(x — vt)], as in the Landau theory, but, in contrary to it, these terms are related to non-damping electrical 
fields. 

6. The spreading of transverse electromagnetic waves in a collisionless plasma is traditionally controlled by the 
Landau rule of passing around the pole and then solving an artificially constructed dispersion equation. But as 
above, in this case the precise (i.e. without Landau additives) dispersion equation for the linearized and Maxwell field 
equations, resulting after substitution therein the travelling wave exp(— iwt + ikx), has no solutions. 

The substitution into the Maxwell field equation of the solution of the kinetic equation, which had been obtained 
with the method of characteristics, leads in non-relativistic approach (after omitting poles terms of the order (v/c) 2 ) 
to the common optical wave equation for a wave in refractive medium with the refracting index n = y/l — lu^/uj 2 ,, 
where is the electron Langmuir frequency, so the solution at real n is represented by the travelling wave 
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E± ~ cxp(— iuj t + ikx) 



(31) 



These results are well known, but it is worth to note however that in the case of medium separating (boundary) 
surfaces there is possible a peculiar non-orthodox solution of the optical equation, as a differential equation with 
ruptured boundary conditions, different from the waves (31), but that one is unstable and questionable to be realized 
experimentally, e.g. in femtosecond spectroscopy [7]. The pole problems do not arise in this consideration at all, 
whatever singularities in the non-relativistic approach are absent (this corresponds to the well known fact, that there 
are large v's in poles, so that the Maxwell distribution function fo(v) there is very small). The pole problem certainly 
arises in the relativistic case and could be solved by applying the above method of two-dimensional transformation to 
the relativistic kinetic equation. 

As for the experimental verification of the Landau theory, first of all we note that there is a very small number of 
such works, which also were carried out by the same scientists [8,9]. Such experiments in fact should be very delicate. 
Certainly, it is very difficult to obtain experimentally a collisionless plasma with Maxwellian electron distribution 
function due to conflicting demands: the Maxwell distribution is just commonly a consequence of electron energy 
exchange in collisions, that is the plasma must be collisional. There is also a possibility in principle of a dependence 
of perturbation spreading in plasma on geometric, spectral and other features of an excitation source. 



Conclusions 

From an unprejudiced consideration of the Landau theory of spreading and damping of plasma waves there is 
emphasized fundamental fact that the dispersion equation, derived in accord with the classical canons, does not have 
whatever solutions, but which must exist at least for the most slowly damping Landau's asymptotic travelling wave, 
be it a case of collisionless or low-collision plasmas. The bypass way, which had been chosen by Landau in 1945 and is 
now cited in all textbooks on plasma physics almost without variations (cf. [4,5]), respectively Landau's rule of passing 
around poles, appears to be unsatisfactory because at tending the imaginary part S of the frequency w to zero resulting 
expressions lose their sense since they do not tend to any definite limits. Thus Landau's analytical continuation is really 
not existent. This is related to the impossibility (e.g. after integration in the dispersion equation by velocity spherical 
angles dipdQ) of expansion in small 8/ojq (that results in asymptotically divergent series, cf. [2]). Landau's theory 
results in a quite sophisticated and mysterious theory of the plasma echo, which leads to the paradox of conflicting 
simultaneous existence of a non-damping mode of plasma echo relative to the Landau field-damping wave. At the 
same time it is possible to construct some reasonable, non-contradictory and mathematically unequivocal theory of 
spreading of plasma perturbations, using a more strict procedure of the two-dimensional Laplace transformation, 
which leads asymptotically to a system of coupled damping and non-damping modes of plasma waves and oscillations, 
but not to an asymptotic single travelling wave with a definite dispersion law, as it should follow from the Landau 
paradigm. 

The rejection of Landau's theory will open many new perspectives in numerous problems of the plasma physics. 

One can explain existing now widely adherence to the Landau theory not only by some natural conservatism and 
piety to Landau, but more by the fact of rather imperceptible, practically unobservable effects arising from this theory 
in all surrounding us realities (laboratory, technical, experimental, observable cosmic, etc.) in spite of the ultimately 
wide its diffusion in abstract sophistications of theorists. 



1 L.D. Landau, J. Phys. (USSR), 2, 25 (1945). 

2 V.N. Soshnikov, V.Ya. Nedzvetsky, Phizika Plasmy, 14, 10, 124S (1988)(in Russian; the journal is translated in English). 

3 I. Alexeff, M. Rader, Int. J. Electronics, 68, 3, 385 (1990). 

4 P.C. Clemmow, J.P.Doughcrty, "Electrodynamics of Particles and Plasmas", 2nd ed. (Addison- Wesley. New York, 1990). 

5 A.F. Alexandrov, L.S. Bogdankevitch, A.A.Rukhadze, ' 'Foundations of Electrodynamics of Plasma", (Vysshaya Shkola, 
Moscow, 1978) (in Russian; there is a prepared English edition). 

6 G. Korn and T.Korn, "Mathematical Handbook for Scientists and Engineers" (McGraw Hill, New York, 1961). 

7 V.N. Soshnikov, Optika i spectroskopiya, 70, 477 (1991) (in Russian; the journal is translated into English). 
8 I.H. Malmberg, OB. Wharton, Phys. Rev. Lett., 17, 175 (1966). 

9 I.H. Malmberg, OB. Wharton, W.E. Drummond, Plasma Physics, 1, 485 (1966). 



7 



